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Abstract. We show how a matrix version of the Buslaev-Faddeev- 
Zakharov trace formulae for a one-dimensional Schrodinger operator 
leads to Lieb-Thirring inequalities with sharp constants 1™ d with y > 
3/2 and arbitrary d > 1 . 



0. Introduction 
Let us consider a Schrodinger operator in L 2 (M d ) 

(0.1) -A + V, 



where V is a real-valued function. In [ |22| ] Lieb and Thirring proved that if 
y > max(0, 1 — d/2), then there exist universal constants L y d satisfying!] 



(0.2) 



tr 



-A + V)I < L Y>C 



Vl + Mxldx, 



In the critical case d > 3 and y = the bound ( P.2[ ) is known as the Cwikel- 
Lieb-Rozenblum (CLR) inequality, see [0, |19|, ^ and also @> 01 • For the 
remaining case d = 1 , y 



1 /2 the estimate ( p.2[ ) has been verified in [|2q], 
see also 0130, On the other hand it is known that ( ]072| ) fails for y = if 
d = 2 and for < y < 1 /2 if d = 1 . 

If V G L Y+ 2 (M d ), then the inequalities Q0.2| ) are accompanied by the 
Weyl type asymptotic formula 



lim 



(0.3) 



1 



tr(-A + aV)I 



a' 



lim 



cl 

~y,d 



2 J 



(l^| 2 + aV)l 



S d xK d 



dxd£, 
(27r) d 



V! +f dx, 
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'Here and below we use the notion 2x_ := |x| — x for the negative part of variables, 
functions, Hermitian matrices or self-adjoint operators. 
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where the so-called classical constant l™ d is defined by 
(0.4) If* = [In)-' ' 



n-\yv, = r(y + 1) H , y >o 



It is interesting to compare the value of the sharp constant L Y d in ( p.2| ) and 
the value of L5, d . In particular, the asymptotic formula ( |031 ) implies that 

(0-5) L^d < L Y>d 

for all d and y whenever ( |0!2| ) holds. Moreover, in [jl|] it has been shown, 
that for a fixed d the ratio L Y d /L^ d is a monotone non-increasing function 
of y. In conjunction with the Buslaev-Faddeev-Zakharov trace formulae 
[0, § one obtains [j^] 



(0.6) L Y , d = I_£ d 

for 

(0.7) d = 1 and y > 3/2 . 

On the other hand one knows that 

T Cl , T 

L Y,d ^ L Y,d 

if d = 1 and 1/2 < y < 3/2 (see [|22[]) or y < 1 and d G N (see EJTTJ] ) . 
Up to now ( PTTj ) was the only case where ( p.6D was known to be true for 



general classes of potentials V G L T+ 2 . Notice, however, that ( J0T6D has 
been proven for various subclasses of potentials. If, for example, O C M d 
is a domain of finite measure, 



V(x) 




as xgQ 

as x G M d \ O ' 



then the equality < JDT5| ) with y = can be identified with the Polya conjec- 
ture on the number of the eigenvalues less than <x for the Dirichlet Laplacian 
in O. It holds true for tiling domains [ |23| ] and has been justified in Jl5| ] for 



certain domains of product structure by using the method of "lifting" with 
respect to the dimension d which is also one of the main ideas of this paper. 
If y > h then ( P.6D is true for arbitrary Cl. This is a simple corollary of 
the Berezin-Lieb inequality (see [0], §5; [ |2"0| ] and also [fT7j]). This approach 
has been extended in [[15]] to the Dirichlet boundary value problems for ma- 
trices of pseudodifferential operators in M d with constant coefficients. The 



Berezin-Lieb inequality was also used in [ [161 ] in order to improve the Lieb 
constant QT9| ] in the CLR inequality for the subclass of Schrodinger oper- 
ators whose potentials are equal to the characteristic functions of sets of 
finite measure. 
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Another example is given in JJp, where the identity ( p.6[ ) with y > 1 and 
d G N has been verified for a class of quadratic potentials. 
We note, that with the exception of (p?7|), the sharp value of L Y> a has been 



recently found in [ 13], where it was proved that for d = 1 and y = 1 /2 



Li/2,1 = 21-1/2! = V2. 

In particular, in higher dimensions d > 2 the sharp values of the constants 
L Yi d have been unknown. 

The main purpose of this paper is to verify ( p.6[ ) for any y > 3/2, d€ N 
and any V G U+7 (M d ). 

In fact, this result is obtained for infinite-dimensional systems of 
Schrodinger equations. Let G be a separable Hilbert space, let 1 g be the 
identity operator on G and consider 

(0.8) -A(g)^ G + V{x) , xeR d , 

in L 2 (IR d , G). Here V(x) is a family of self-adjoint non-positive operators 
in G, such that tr V G L Y+ ~ (IR d ). Then we prove that 



(0.9) tr (-A®1 G + V(x))I<lf d 



trVl + Mx)dx 



for all y > 3/2 and d > 1 . The inequality ( p.9[ ) can be extended to magnetic 
Schrodinger operators and we apply it to the Pauli operator. 

We shall first deduce (£9} for d = 1, y = 3/2 and G = C n from 
the appropriate trace formula (|1.61|) for a finite system of one-dimensional 



Schrodinger operators. In the scalar case these trace identities are known as 
Buslaev-Faddeev-Zakharov formulae []5J,[8|]. The matrix case can be handled 
in a similar way as in the scalar case (see [§]) However, we give rather 
complete proofs of the corresponding statements in section 1, since we were 
unable to find the necessary formula ( [1.61D in the numerous papers devoted 
to this subject. 

Note that we discuss trace formulae only as a technical tool in order to 
establish bounds on the negative spectrum. We therefore develop the theory 
of trace identities only as far as it is necessary for our own purpose. 

In section 2 we extend the results of section 1 to the Schrodinger operator 
in L 2 (M 1 , G). Applying a "lifting" argument with respect to dimension as 



used in [Bp and [15], we obtain in section 3 the main results of this paper. 



Finally we would like to notice that the combination of the results of this 
paper and the equality 1-1/2,1 = 1/2 discovered in QT~3| ] has lead to new 
bounds on the Lieb-Thirring constants in [|T2p which improve the corre- 
sponding bound obtained in [Q] and [ |2"T| ] . 
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1. Trace formulae for elliptic systems 

1.1. Jost Functions. Let and 1 be the zero and the identity operator on 
C n . We consider the system of ordinary differential equations 

(1.1) -Q^®l)yM+V(x)y(x) =k 2 -y(x), xGl, 

where V is a compactly supported, smooth (not necessary sign definite) 
Hermitian matrix-valued function. Define 

x m i n := min supp V and x max := max supp V . 

Then for any k £ C\{0} there exist unique n x n matrix-solutions F(x, k) 
and G (x, k) of the equations 

(1.2) -F^(x,k)+VF(x,k)=k 2 F(x ) k), 

(1.3) -G^(x,k) +VG(x,k) =k 2 G(x,k) , 
satisfying 

(1.4) F(x,k) = e ikx 1 as x > x max , 

(1.5) G(x,k) = e~ ikx 1 as x < x min . 

If k £ C \ {0}, then the pairs of matrices F(x, k), F(x, — k) and G(x, k), 
G(x, — k) form full systems of independent solutions of fll.ip. Hence the 
matrix F(x, k) can be expressed as a linear combination of G(x, k) and 
G(x,-k) 

(1.6) F(x,k) = G(x,k)B(k) + G(x,-k)A(k) 
and vice versa: 

(1.7) G(x,k) = F(x,k)|3(k) + F(x,-k)cx(k) . 

1 .2. Basic properties of the matrices A(k) , B (k) , <x(k) and (3 (k) for real 

k. Throughout this subsection we assume that k £ R\{0}. Consider the 
Wronskian type matrix function 

VV| [F, G](x,k) = G*(x,k)F x (x,k) - (G x (x,k))*F(x,k) . 



Then by (jlj) and (jL3D for k £ R we find that 

A Wl [F, G] (x, k) = G*(x, k)F x (x, k) - (G£(x, k))*F(x, k) = . 
dx 

Note that for x < x min by ( |TT6| ) we have 
Wi [F, G](x,k) = [G*(x,k)G x (x,k) - (G' x (x,k))*G(x,k)] B(k) 

+ [G*(x, k)G x (x, -k) - (G'Jx, k))*G(x, -k)] A(k) 
= -2ikB(k) , 
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while for x > x max by ( |T77| ) we find 

Wi[F, G](x,k) = B*(k) [F*(x,k)F x (x,k) - (F x (x, k))*F(x, k)] 

+ <x*(k) [F*(x,-k)F x (x,k) - (F x (x,-k))*F(x,k)] 
= 2ik(3*(k) . 
This allows us to conclude that 

(1.8) 6*(k) =-B(k) . 
Similarly, for the matrix- valued function 

W 2 [F, G](x,k) = G*(x,k)F x (x,-k) - (G x (x,k))*F(x,-k) 
we have ^W 2 [F, G] (x, k) = and 

W 2 [F, G](x, k) = -2ikA(-k) as x < x min , 

W 2 [F, G](x, k) = — 2ikcx*(k) as x > x max . 

Thus, 

(1.9) A(-k) = <x*(k). 

Inserting ( |1.6| ) into ( |1 .7[ ) and making use of (|TTE|), ( |1.9| ) we obtain 

(1.10) G(x,k) = G(x,k) [B(k)|3(k) + A(-k)cx(k)] 

+ G(x,-k)[A(k)(3(k)+B(-k)a(k)] , 

and thus 

(1.11) A(-k)A*(-k) -B(k)B*(k) = 1 , 

(1.12) B(-k)A*(-k) -A(k)B*(k) =0. 
In particular, this implies 

(1.13) |detA(k)| 2 = det A(k) det A*(k) =det(l + B(-k)B*(-k)) > 1 
for all k G R\{0}. 

1.3. Associated Volterra equations and auxiliary estimates. Next we 
derive estimates for the fundamental solutions of ( |1.1| ) for Im k > 0. Note 
first that the matrices F(x, k) and G(x, k) are solutions of the integral equa- 
tions 



(1.14) F(x,k) = e ikx 1 



(1.15) G(x,k) = e- ikx 1 + 
Put 

H(x,k) = e~ i,oc F(x, k) - 1 



k _1 sin k(x - t)V(t)F(t, k) dt , 
k" 1 sin k(x - 1) V(t) G (t, k) dt 
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Obviously, this matrix-valued function satisfies 



(1.16) 
and 

(1.17) 
where 
(1.18) 

Note that 
(1.19) 



H(x, k) = for x > x n 



H(x,k) 



K(x,t,k) dt + 



K(x,t,k) 



,2ik(t-x) 



2ik 



K(x,t,k)H(t,k) dt, 



-V(t) . 



|K(x,t,k)|| < Ci(V,n)/(1 +|k|) 



for all k with Im k > and all k with x m i n < x < t. Here and below | 
denotes the norm of a matrix on C n . 

Solving the Volterra equation Q1.17D we obtain the convergent series 



H(x,k) = ^ 



m=l 



] J K(x l _ 1 , x b k) dx-i ■ ■ ■ dx r 



X<X] <-<x m 



1-1 



From ( |1.19[ ) we see that |H(x, k)| < C2CV) for all x m ; n < x < x max . Insert- 
ing this estimate back into ( J1.17| ), we conclude that the inequality 



(1.20) 



|H(x,k)|| < C 3 (V,n)(1 +|k| 



holds for all x with x min < x < x max and all k with Im k > 0. 

Remark 1.1. If we assume that Im k > and |k| > 1 , then ( ]1.19D and 

therefore ( |1 .20[ ) holds true for all x£K. 

It is not difficult to observe, that H(x, k) defined by Q1.17D is smooth in 
(x, k) G R x {k G C : Im k > 0} . 
In particular, if we differentiate ( |1.17| ) with respect to k we find that 



T=H(x,k) 
ok 



K(x,t,k)^H(t,k) dt 
3k 



Since 3H(x, k) /3k satisfies a homogeneous Volterra integral equation with 
the kernel ( |1.18D , we obtain 3H(x, k)/3k = 0, and thus all the entries of 
the matrix H(x, k) are analytic in k for Im k > 0. 
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1.4. Further Estimates on A(k) andB(k). If we rewrite ( |1 . 14p as follows 



(1.21) F(x,k) = e 



+ 



ikx 



1 

2ik 



V(t) dt 



1 

2ik 



V(t)H(t,k) dt 



,— ikx 



2ik 



3 2ikt V(t) dt + 



e 2ikt V(t)H(t,k) dt 



then the expressions in the brackets in the r.h.s. do not depend on x for 



x < x m i n . Comparing ( |1.21[ ) with ( |L6] ) we see that 



(1.22) A(k)=l 



(1.23) 



1 

21k, 

r+oo 



V(t) dt 



1 

2ik. 

e 2ikt V(t) dt + -L 
2ik 



V(t)H(t,k) dt, 



oo 
r+oo 



e 2ikt V(t)H(t,k) dt. 



For sufficiently large |k| > C the smoothness of V and Q1.20| ) imply 



(1-24) 



A(k) 



1+ 2Tk 



V(t)dt 



< C 4 (V,u) |k 



-2 



Im k > , 



(1.25) ||B(k)|| < C 5 (V,n)|kr 2 , kGl. 
In subsection 1.6 we shall see that ( |1 .25[ ) can be improved so that 

(1.26) B(k)=0(|kr m ) for all meN as k -> ±oo . 

1 .5. The matrix A (k) for Ira k > 0. First note that all entries of the matrix 
A(k) are analytic in k for Im k > and continuous for Im k > 0, k ^ 0. 
This follows from Q1.22P and the analyticity of H(x, k) . Fixing a sufficiently 
small e > and by using ( |1 .22[ ) and ( |1.20[ ) we obtain 

(1.27) ||A(k)|| < Cglkr 1 as |k| < e, lmk>0. 

Moreover, all the entries of A(k) and thus the function det A(k) are analytic 
for Im k > and continuous for Im k > 0, k 7^ 0. Near the point k = we 
find 



(1.28) 



|detA(k)| < C 7 |k|" 



as 



|k| < e, lmk>0. 



Next let us describe the connection between the function det A(k) and 
the spectral properties of the self-adjoint problem ( |1.1| ) on L 2 (M, C n ). Our 
assumptions on the matrix potential V imply, that the operator on the l.h.s. 
of fll.lD has a discrete negative spectrum, which consists of finitely many 
negative eigenvalues \\ = (ixt) 2 , k\ > of finite multiplicities mi. Ob- 



viously a solution y(x) of fll.lD with k = \h\ belongs to L 2 ' 
only if 

G(x, ixje^ as 



L,C n ), if and 



y(x) = F(x, ix\)e Y y as x > x n 



8 A. LAPTEV AND T. WEIDL 

for some non-trivial vectors e G , ep 6 C n . Linear independent so- 
lutions y i , . . . , y m i define linear independent vectors , . . . , e^ m and 
, . . . , , respectively. In view of (jL^) we conclude that 

(1.29) dimker A(ix;) = mi . 

If we select an orthonormal basis in C n , such that the first mi elements 
belong to ker A(ixi), we find that the first mi rows of A(k) vanish as k — > 
Since det A(k) does not depend on the choice of the orthonormal basis 
and all entries of A(k) are analytic, the function det A(k) has a zero of the 
order 

(1.30) m' l >m l 

at k = \K\, H\ > 0. Moreover, if A = k 2 , Im k > is not an eigenvalue of 
the problem flOD, then det A(k) ^ 0. 

In the remaining part of this subsection we prove that 

(1.31) m^ = mt. 

Let g (x, y , k) be the Green function of the problem Ql.l|). If k 2 < 0, Im k > 
0, and det A(k) ^ it can be written as 



g(x,t),k) 



G(x,k)Z (y,k) as y > x 
-F(x, k)Z + {y, k) as y < x 



Here Z + (y , k) and Z [y , k) are n x n-matrices, which are chosen such that 
lim g(x,y;k) = lim g(x,y;k) , 

x=y— x=y+0 

lim g^(x,y;k)= lim g' x (x,y;k) + 1 . 

x=y— x=y+0 

These equations turn into 
(1.32) 

wt».«©a-(?). ®)- 

Since ^ det W = 0, the determinant of W is a constant with respect to y . 
If y with y < x m j n , Im k > 0, then in view of ( |1.6D and ( JT3| ) we have 

(133) Wfu k) = ( e-^B(k) + e^A(k) \ 

K } VV[ y> K > ^-tice-ikvi -ike- iky B(k) +ike iky A(k) y l ' 

Hence 

detW= (2ik) n detA(k) 
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and W is invertible if and only if det A(k) ^ 0. From fll.33| ) we see then, 
that for y < x min the entries of 



(1.34) W-\y,k) 
satisfy 



Xn(y,k) X 12 (y,k) 
X 21 (y,k) X 22 (y,k) 



e~ iky X 21 - ike~ iky X 22 = , 
e~ iku (X 21 - ikX 22 )B(k) + e iky (X 21 + ikX 22 )A(k) = 1 . 

This gives X 21 (y , k) = ikX 22 (ij , k), and thus 

X 22 (-y,k) = (2ik)- 1 e- iky A- 1 (k). 

In view of ( |T32| ) and (jO|) we obtain Z+(y,k) = X 22 (-y,k) and finally 
conclude that 

(1.35) 

g(x,y,k) = -(2lk) _1 A _1 (k)e lk(x_v] as y < x min < x max < x . 

If k is in a sufficiently small neighbourhood of ixj., the Green function 
g (x, y , k) can be written as 



g(x,y,k) = — =H — — — — — - + g 1 (x,y,k) . 
(k- xx0(k + ixO 

Here gi(x,y,k) is locally bounded and {ihrl^i f° rms an orthonormal 
eigenbasis corresponding to the eigenvalue Ai = — k\. Hence, 

detX 22 (-y,k) = (2ik)- n e^ lnku detA- 1 (lc) 

= (-l) n e- inkx detg(x,y,k) = O (|k - ix^™ 1 ) 

as k — > ixi. This implies that det A(k) has a zero of the order 

at k = \K\, Finally, the last inequality and ( |T3| ) imply ( |1.31| ). 
1.6. The matrix function T(x, k). Consider the matrix function 



(1.36) T(x,k) = 1 +H(x,k) = 1 + 



K(x,t,k)T(t,k) dt 



According to subsection 1.3 the matrix- valued function T(x, k) is smooth 
and uniformly bounded for 

(x,k) G R x {k G C : Im k > and |k| > 1} . 
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Obviously T(x, k) = 
using ( |1.18| ) we obtain 
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1 for x > x max . Integrating by parts in ( |1.36D and 



(1.37) 



dx* 



T(x,k) 



,2ik(t-x) 



— (V(t)T(t,lc))dt 

for all I e N. Since supp V C [x min , x max ] we find 

(1.38) d l T(x,k)/dx l =0 as x max < x , 

(1.39) ||d l T(x,k)/dx l || < C 8 as x min < x < x max , 

(1.40) ||d l T(x,k)/dx l || <C 9 e 2(x - Xm '" )Imk as x < x min , 

for all k with Im k > and |k| > 1 . The constants Cs and C9 depend only 
upon V, n and I. If we integrate the r.h.s. of ( |1.37[ ) by parts, then fll.39j ) 
and ( |1 .40[ ) imply 

Cio 



(1.41) ||d l T(x,k)/dx l || < 



1+|k| 
Ci 



as x m ; n ^ x ^ x max , 



(1.42) ||d l T(x,k)/dx l || < — il^ e 2(x - x ™" )Imk as x < x min , 

11 11 1+|k| 

for all k with Im k > and |k| > 1. The constants Cio and Cn depend 
only upon V, n and I. 

In a similar way integrating by parts in ( |1.37| ), we obtain the asymptotical 
decompositions 



r°o a\— 1 

2ik(t-x^ a 



— (V(t)T(t,k))dt 



\T+~\ AT+l-2 



(1.43) 



— (2ik) r dx r+i 
+ (-l) q+1 
1 



(V(x)T(x,k)) 



'oo _2ik(t— x) jq+l— 1 

-— . F— ri (V(t)T(t,k)) dt 



f 1AT+1 AT+l-2 



L 

r=1 



(2ik) r dx 



r+l-2 



(V(x)T(x,k))+R q>l (x,k) 



as 



as |k| > 1 , Im k > 0. Here 

(1.44) R q ,i(x,k) =0 

(1.45) ||R q>l (x,k)|| <C U (1 +|k|)- q 

(1.46) ||R q ,i(x,k)||< ^ e 2 ^^ k as 
The constants Cu and C13 depend upon V, n, I and q. 



"max — " 



aS X mm ^ X ^ "max > 



X X1-1 
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Since d l H/dx l = d l T/dx l for all I G N, integration by parts in ( |1.23| ) 
and the inequalities ( |T38| ), ( |L4T| ) and ( |L42| ) give ( |L26| ). 



1 .7. The matrix function cr(x, k). By using Q1.16D , ( |1 .20[ ) and Remark 1 . 1 
for sufficiently large |k|, Im k > 0, the matrix T(x,k) = 1 + H(x,k) is 
invertible for all x 6 R and 

(1.47) ||T _1 (x,k)|| < C 14 for all xel, |k|>C 15 , Im k > , 

with sufficiently large constants C14 = Ci4(V,n) and C15 = C-\s{V, n). 
Hence, for sufficiently large |k| with Im k > the matrix function 



(1.48) 



cr(x, k) 



dx 



T(x,k) 



T _1 (x,k) 



is well defined for all xei Liouville's formula 



— (lndetT(x,k)) = tr 
dx 



i-T(x,k) 
dx 



T-^x.k)! 



implies 



dx 



(lndet e~ ikx F(x, k)) = tr o(x, k) 



Since e tkx F(x, k) = 1 as x > x max and 



-,— ikxi 

[K, kj — 1 as A. ^ A. max 

e- ikx F(x,k) = e~ 2ikx B(k) + A(k) = A(k) + o(l) 
as x ^ — 00, Im k > e > 0, we finally conclude that 



(1.49) 



lndetA(k) = 
|k| > C 15 , 



tr cr(x, k) dx , 

—00 

Im k > e > . 



Remark 1.2. Formula Q1.49D is a matrix version of the corresponding well- 
known identity for scalar Schrodinger operators (see e.g. §3 in [§]). 



1.8. The asymptotical decomposition of cr(x, k). Next we shall develop 
cr(x, k) into an asymptotical series with respect to the inverse powers of k. 
For the sake of future references we compute the first three terms, although 
we only need the second one in this paper. 

If we apply ( |1.43| ) with q = 2, 1 = 1 we find that 



(1.50) 



Q2 = Rz.iF- 1 , 
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while fll .43| ) with q = 4, 1 = 1 gives 
(1.51) a 1 



d 2 V dV d 2 T , 
(2ik) 3 dx~ ff+ dx 2 

"(2ib{^ + Vff } + 2lk V+R4 ' ,T "- 



Inserting ( |1.50| ) into ( |1.51[ ) we obtain 



(1.52) 



1 dV 



1 



2ik (2ik) 2 dx (2ik) 3 



V 2 



d 2 V 
dx 2 



+ Q 4 - 



Finally, if we insert in a similar way ( |1.52| ) and ( |1 .43[ ) with I = 2, q = 3 as 
well as I = 3, q = 2 into ( |1 .43[ ) with I = 1 and q = 6, we arrive at 

,dV , f d 2 V ,1 

(1.53) a = (2ik)- 1 V - (2ik)- 2 — + (2ik)- 3 1 — -V 2 l 



- (2ik)- 4 
+ (2ik)- 5 



d 3 V 2 dV 2 



dx 3 

d 4 V 
dx4 



dx 

d 2 V 2 
' dx 2 



S +2V3 



Q* 



As well as R q> t the terms Qz, Q4 and satisfy the inequalities of the 
type ( |1.44[ ) - ( J1.46D with q = 2, q = 4, and q = 6, respectively. Then we 
conclude that 

p+00 

trQ q (x,k)dx = 0(|k|- q ), q =2,4,6, 



as |k| — > 00 with Im k > e > and thus, 
(1.54) 

1 r+°° "] 

2ikL' rVdX -(2ikF 
1 r+oo 



tr cr(x, k)dx 



r+oo 



tr V 2 dx 



as Ikl 



(2ik) 5 J_ c 
00 with Im k > e > 0. 



2tr V 3 + tr ( — 

dx 



dx + 0(|kj 



-61 



1.9. The dispersion formula. Let 

{AiK N =1 ={(tx l ) 2 }^ 1 , x t >0, 

be the finite set of the negative eigenvalues of (|1.1|). Each eigenvalue occurs 
in this set only once. Let mi be the order of zero of det A(k) at the point 
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k = ixj., which by section 1.5 equals the multiplicity of the corresponding 
eigenvalue. Then the arguments in section 1 .5 imply that the function 



(1.55) 



M(k) =hJ detA(k)]^J 



i-i 



k + XKx 

k — \k\ 



is analytic for Im k > and continuous up to the boundary except k = 0, 
where it has at most a logarithmic singularity. Moreover, the inequality 
(|1.24p gives 



|M(k)| < C 2 (V)|k| 



-i 



for all sufficiently large |k| > C, Im k > 0. Hence, by applying Cauchy's 
formula for large semi-circles in the upper half-plane we obtain 



+oo 



M(z) dz 



(2m)M(k) , 



M(z) dz 







for arbitrary k with Im k > 0. This implies 

00 Re Mfz) 



(1.56) 



M(k) 



1 

ni 



dz, 



which by ( |1.55[ ) is equivalent to 



(1.57) lndetA(k) 



7tl 



N 



In | det A(z)| dz v— k — ix^ 

h > m t ln- 

z — k z — 



k + ix\ 



for all k with Im k > 0. 



1.10. Trace formulae for elliptic systems. Note that 



T— Y-XXi 2 t— I v- 3 

) itu In t—^ — = > m, X Hx - ——- ) m l x 1 



i=i 

N 



(1.58) 



i=i 



as |k| — > co, Im k > e > 0. On the other hand from Q1.13D and ( |1.26D we 
have 



ln|detA(z)| = 2- 1 ln|det(l + B(-z)B*(-z))| = 0(|zr m ), zgK, 
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as |z| — ) oo, for all m G N. Hence, the integral in ( |1 .57[ ) permits the 
asymptotical decomposition 



(1.59) 



lnldetAfzll dz 



m 



j=0 

P+OO 



z' ln|det A(z)| dz 



as |k| — > oo, Im k > e > 0. 

Combining ( |1.58| ), fll.59| ) with m = 5 and ( |1.54| ) we obtain 

N 



(1.60) 
(1.61) 
(1.62) 



3 

T6 



tr Vdx 



tr V 2 dx 



i° -r 



2tt 
3I 2 



N 



J5_ 

32 



tr V 3 dx + 



64 



to- 



dy 
dx 



dx = 



2n + Z m ^ 3 ' 

1=1 

N 



5U 
2ti 



i=i 



Finally we remark, that in view of ( |1.13| ) 
(1.63) 

for all even, non-negative integers j 



Ij > 



2. Sharp Lieb-Thirring inequalities for second order 
one-dimensional schrodinger type systems 

2.1. A Lieb-Thirring estimate for finite systems. Let us first consider the 
operator on the l.h.s. of ( |1.1| ) in L 2 (M, C n ) for some smooth, compactly sup- 
ported Hermitian matrix potential V. Preserving the notation of the previous 
section the bounds ( |1.61[ ) and ( |1.63| ) imply 

a2 x 3 / 2 



(2.1) 



tr 



dx 2 



i +v(x; 



Y_ m i x i 



< 



16 



trV 2 (x) dx. 



By continuity (£Oj) extends to all Hermitian matrix potentials, for which 
tr V 2 is integrable. Finally, a standard variational argument allows one to 
replace V by its negative part V_: 



(2.2) 



tr 



dx 2 



i + v(x; 



3/2 



~ 16 



tr V 2 (x) dx. 



The constant in the r.h.s. of this inequality is sharp and coincides with the 
classical constant Ly 2 r In particular, this constant does not depend on the 
internal dimension n of the system. 
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2.2. Operator-valued differential equations. Let G be a separable 
Hilbert space with the scalar product < •, • >g and the norm || • ||g- Let 
H 1 (R, G) and H 2 (R, G) be the Sobolev spaces of all functions 



u 



for which the respective norms 



MIhi 



U 



Iu'IIg + IHg) dx 
|u"|| 2 g + H| 2 g) dx 



are finite. Finally, let 1 g be the identity operator on G. Then the operator 
~~ cE? ® ^ g defined on H 2 (R, G) is self-adjoint in L 2 (R, G). It corresponds 
to the closed quadratic form 



h[u, u] 



|u || G dx 



with the form domain H 1 (R, G). 

Let B and % respectively be the spaces of all bounded and compact lin- 
ear operators on G. Let || ■ || 2 denote the corresponding operator norm. 
Consider an operator- valued function 



W(-): 



for which W(x) = (W(x))\ x G Rand ||W(-)||s e L p ( 
Denote 



1 < p < oo. 



w[u, u] 



(W(x)u(x),u(x)) G dx. 



This form is well-defined on H 1 (R, G) and 



(2.3) 



|w[u,u]| < C 16 



||W(x)||^dx ||u 



Vp 



H 1 



The constant Ci6 does not depend upon W or u. Moreover, for all e > 
there exists a finite constant C-\j{e, W), such that 



(2.4) 



|w[u,u]| < eh[u, u] + Ci7(e, W) 



\u\\r dx . 



Both ( |2.3[ ) and ( |2.4p follow immediately from the corresponding inequali- 
ties which hold in the scalar case. Hence, the quadratic form 



h[u, u] + w[u, u] 
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is semi-bounded from below and closed on H 1 (R, G). It induces a self- 
adjoint semi-bounded operator 

(2.5) Q = —^- (g) i G + wix) 

on L 2 (R,G). 

If in addition W(x) £ % for a.e. x £ R, then the form w[-, •] is relative 
compact with respect to the metric on H^R, G). In order to prove this 
fact we introduce the orthogonal projections Pm on the linear span of the 
first M elements of some fixed orthonormal basis in G. As a consequence, 
the Birman-Schwinger principle implies, that the negative spectrum of the 
operator Q is discrete and might accumulate only to zero. In other words, 
the operator Q is compact on L 2 (R, G). 

2.3. A Lieb-Thirring estimate for operator- valued differential equa- 
tions. We shall prove the following Theorem: 

Theorem 2.1. Let W[x) be self-adjoint Hilbert-Schmidt operators on G 
for a.e. x £ R and let tr W 2 (-) £ L 1 (R, G). Then we have 



3/2 

(2.6) tr + <lf /21 



tr W 2 dx . 



where according to ( |074t ) it holds Ly 2 1 =3/16. 

Proof. Assume that ( ^o] ) fails. Then there exists a non-positive operator 
family W satisfying tr W 2 (-) £ lJ(R) and some sufficiently small e > 0, 
such that 



(2.7) tr X 3 / 2 (Q)>^ 



foo 



tr W z dx . 

— oo 



Here 



Xe (Q) = -E ( _ 00i _ e) (Q)Q ) 



with E(_ 00 _ e )(Q) being the spectral projection of Q onto the interval 
(— oo, — e). Since Q is compact, the operator E(_oo e )(Q) is of a finite 
rank n(e). 

Fix some orthonormal basis in G and let Pm be the projection on the 
linear span of its first M. elements. Consider the auxiliary operators 

Q(M, e) = En.o.-ejfQKI (x) <g> P M )Q0 M <g> Pm)E ( -oo,- £ )(Q) • 

Obviously we have rank Q(M, e) < n(e) for all M. Since 1 (x) ® Pm 
turns to the identity operator on L 2 (R, G) in the strong operator topology 
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as M — > oo, then the operators Q(M, e) converge to Xe(Q) i n me L 2 (]RG' 
operator norm, as M — > oo and 



tr (Q(M, e)) 3 / 2 -> trx e (Q) as M -> oo . 



Thus, 
(2.8) 



tr(Q(M,e))L /2 >^ 



f+OO 



tr W 2 dx 



for some sufficiently large M. On the other hand, a standard variational 
argument implies 

tr (Q(M, e)) 3/2 < tr ((1 (x) ® P M )Q0 W ® P M )) 3/2 • 



Observe that the expression on the r.h.s. is nothing else but the Riesz mean 
of the order y = 3/2 of the negative eigenvalues of the M x M-system (JTTT|) 
with V(x) = P M W(x) Pm- Thus, from (TO) we obtain 



tr(Q(M,e))i /2 < A 



tr V 2 (x) dx < — 
1 6 



trW 2 (x) dx, 



which contradicts Q2.8|). This completes the proof. 



2.4. Lieb-Thirring estimates for Riesz means of negative eigenvalues of 
the order y > 3/2 . We shall now suppose, that the non-positive operator 
family W(x) satisfies 



(2.9) 



trW_ 



fx) e L 



for some y > 3/2 . 



Let dEf-oo AjlQ) be the spectral measure of the operator Q. Repeating the 
arguments of Aizenman and Lieb [JT]], we find 



B(y-^J)trQl = tr 



dE 



(-oo,A) 



t y -2tr (Q + t 



(q: 



,3/2 



t Y ^(t + A) 3/2 dt 



dt 



< 



o 

3 r°° 

16 Jo 



r+oo 



dtt^2 



tr (W(x) + 1) 2 dx , 
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where B(x,ij) = j^w^ is the Beta function. Let — (Xj(x) < be the 
negative eigenvalues of W(x). Then 



dtt^2 



tr (W{x)+t) 2 _dx 



p+oo 



dx 



B(y-2- 3 



dx y~ \i 



B(y-2-3 



r+oo 



tr W_ 



fx) dx . 



From ( P.4| ) we obtain 



r( y + r 



-T.i 2 7tV2r(y + §; 

and this implies 



£ r(y+1)r(3) _ 3 B(y-|,3) 

i6'r( Y + |)r(f) 16'b(y-|,|)' 



Theorem 2.2. Let the non-positive operator family W(x) satisfy ( |2.9| ). 



(2.10) 



tr 



dx 2 



1g + W(x) <I« 



tr Wl + Mx) dx. 



It remains to note, that the constant L~-, in ( |2.10| ) is approached for po- 
tentials aW as cx — > +oo. 

3. LlEB-THIRRING ESTIMATES WITH SHARP CONSTANTS FOR 
SCHRODINGER OPERATORS IN HIGHER DIMENSIONS 

3.1. Lieb-Thirring estimates for Schrodinger operators. Let G be a 

separable Hilbert space. We consider the operator 

(3.1) -A®1 G + V(x) 
in L 2 (M d , G). If the family 

V(-) :R d ^ "B 
of bounded self-adjoint operators on G satisfies 

(3.2) ||V(-)||s € L p (R d ), max{1 , d/2} < p < oo , 
then the quadratic form 



v[u,u] 



(V(x)u(x), u(x)) G dx 
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is zero-bounded with respect to 



h[u, u] 



d 

L 



3u 



3xi 



dx . 



This immediately follows from the corresponding scalar result and the argu- 
ments given when proving the inequalities (|23|), (^3|). Hence the quadratic 
form h[-, ■] + v[-, •] is semi-bounded from below, closed on the Sobolev 
space H^IR* 1 , G) and thus generates the operator ( |3.1| ). As in subsection 
3.2 one can show, that if in addition to ( p.2| ) we have V(x) G % for a.e. 
x G IR d , then the negative spectrum of the operator (^jj) is discrete. 
The main result of this paper is 

Theorem 3.1. Assume that V(x) < Ofor a.e. x G M d and that tr V^ +Y (-) 
is integrable for some y > 3/2. Then 



(3.3) 



tr 



-A®1 6 + V(%))I<I* d 



trVi +Y (x) dx. 



Proof. We use the induction arguments with respect to d. For d = 1 , 
y > 3/2 the bound ( |3.3D is identical to ( |2.10| ). Assume that we have ( |3~3| ) 
for d — 1 and all y > 3/2. Consider the operator ( p.l[ ) in the (external) 
dimension d. We rewrite the quadratic form h[u, u] + v[u,u] for u G 
H 1 (E d ,G) as follows 



r+oo 



H[u, u] + v[u,u] 
h(x d ) [u, u] 

w(x d )[u,u] 



h(x d ) [u, u] dx d + 



w(x d ) [u, u] dx d , 



3u 



3x d 

d-l 

L 

j=1 



dX! • ■ ■ X d _T , 



3u 



3xi 



+ (V(x)u,u) ( 



dxi • ■ ■ X d _! 



The form w(x d ) is closed on H 1 (M d 1 , G) for a.e. x d G R and it induces 
the self-adjoint operator 



W(x d ) =~Y_ 9^2 ® lG+V(x 1 ,...,x d _ 1 ;x d ) 

on L 2 (R d_1 ,G). The negative spectrum of this (d — 1 ) -dimensional 
Schrodinger system is discrete, hence W_(x d ) is compact on L 2 (IR d ~ 1 , G) 
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and according to our induction hypothesis tr W_ 2 (x d J satisfies the in- 
equality 

(3.4) 



•y_| 

tr V 2 (xi , . . . , x d _-,; x d ) dx! • • • dx d _! 



for a.e. x d £ R. For V £ U+7(R d - 1 ), the function tr Wl +2 {■) is inte- 
grable. 

Let w_(x d )[-,-] be the quadratic form corresponding to the opera- 
tor W_(x d ) on H = L 2 (R d- \G). Then we have w(x d )[u,u] > 
— w_(x d )[u, u] and 



(3.5) ri[u, u] + v[u, u] > 



3u 



3x r 



(W_(x d )u,u) H 



H 



dx d 



for all u £ H 1 (R d , G). According to section 2.2 the form on the r.h.s. 
of ( |3.5p can be closed to H 1 (R, H) and induces the self-adjoint operator 

® 1 H - W_(x d ) 



dx 2 



on L 2 (R, H). Then (g3J) implies 

( d 2 V 
(3.6) tr (-A ® 1 G + V)I < tr (--^ ® 1 H - W_(x d ) J . 

We can now apply ( [2.10D to the r.h.s. of ( ]3T6| ) and in view of Q3.4| ) we find 
d 2 



tr 



dx 2 



1 H -W_(x d )) 



tr Wl +2 (x d ) dx d 



< T Cl T Cl 



trVZ +2 (x)dx 



The calculation 



L^iL 01 



r(y + i) 



y+T ,d-1 



27tir(y+i + l) 2 d - 1 7r ii T i r(y + ^ + ^ 1 + l 



r( y + i 



= L 



cl 



2 d 7tfr(y + | + 1 

completes the proof. 

For the special case G = C we obtain the Lieb-Thirring bounds for scalar 
Schrodinger operators with the (sharp) classical constant L y d = Ur, d for 
y > 3/2 in all dimensions d. 



LIEB-THIRRING INEQUALITIES 



21 



3.2. Lieb-Thirring estimates for magnetic operators. Following a re- 
mark by B. Helffer [ ]T0| ] we demonstrate, how Theorem [3J] can be extended 
to Schrodinger operators with magnetic fields. Let 

a(x) = (a 1 (x),...,a d (x)) t , d>2, 

be a magnetic vector potential with real- valued entries G L 2 oc (IR d ). Put 

H(a) = (iV + a(x)) 2 ®1 G . 

Its form domain d[h(a)] consists of the closure of all smooth compactly 
supported functions with respect to H(a)[-, ■] + 1| • \\12ma G 



(cf. [E5]]), where 



h(a)[u,u] = Y_ 



k=1 " 



3xi 



a k u 



dx , 



Let the operator family V and the corresponding form v be defined as in 
the previous subsection. If ( JO ) is satisfied, then one can apply Kato's 



inequality [ Jl4| , |25| ], and find that the form 

(3.7) q(a)[u,u] = h.(a)[u,u] +v[u,u] 

is closed on d[q(a)] = d[H(a)] and induces the self-adjoint operator 

(3.8) Q(a) =H(a)+V(x) 

on L 2 (R d , G). Finally, by applying Kato's inequality to the higher- 
dimensional analog of (pT3|) we see, that V(x) G % for a.e. x G M d in 



conjunction with ( J3T2| ) implies that Q ( a) has discrete negative spectrum. 

Theorem 3.2. Assume that a G L 2 oc (IR d ) is a real vector field, and that 
the non-positive operator family V(x) satisfies tr V2 +y g L 1 (M d ) for some 
y > 3/2. Then 



(3.9) 



tr (H(o)+V(x))I<e d 



trVi +Y dx. 



Proof. In the dimension d = 1 , any magnetic field can be removed by gauge 
transformation. Thus ( |2.10[ ) can serve to initiate the induction procedure. 

Assume now that ( |3.9[ ) is known for all y > 3/2 for the dimension d — 1 
and consider the operator H(a) in the dimension d. Put 



W(x d ) = 



d-1 

L 

.n=1 



3x T 



+ a T 



V(x) 
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> 











2 


dx + 


R d 


(i 9 

V 9 *d 




u 


G 










2 


dx — 


R d 


(i 9 
V 9x d 


+ a d ^ 


u 


G 



(W(x d )u,u) H dx d 
(W_(x d )u,u) H dx d , 



where for fixed x d G R the operator W_(x d ) is the negative part of W(x d ) 
on H = L 2 (M d ~ 1 , G). We now choose a gauge in which a d vanishes. 
Namely, put 



4>(X!,...,X d ) 



a d (x 1 ,...,x d _ 1) T)dT 



andu(x) = e i<i)M u{x) for all u G d[q(a)]. Then 
(3.10) 







3u 


2 


> 






dx — 




R d 


3x d 


G 



(W(x d )u,u) H dx d , ued[q(a)], 



where 



W(x d ) = e i * (x '' Xd) W_(x d )e- i * (x ' 



[Xt, . . . ,X d _-jJ , 



acts on H for any fixed x d G R. Closing the form on the r.h.s. of ( |3.1C| ) we 
see that 



(3.11) 



tr (H(a)+V(x))I<tr 



d 2 ~ V 

lH-W(x d )J , 



dx 2 



where the operator on the r.h.s. acts in L 2 (R, H). From our induction hy- 
pothesis we have 



trVl +2 (x';x d ) dx'. 
Hence ( |2.1U| ) can be applied to estimate the r.h.s. of ( p.ll| ) and we complete 



trW^(x d ) =tr Wr j lx d ) < L^, d _ 1 



the proof of (pT9|) in the same manner as in the proof of Theorem [3T 



3.3. Lieb-Thirring estimates for the Pauli operator. As an application 
of Theorem |3]2| we deduce a Lieb-Thirring type bound for the Pauli opera- 
tor. Preserving the notations of the previous subsection we put d = 3 and 
G = C 2 . Let a(x) = (ai(x), a2(x), a^x)) 1 be a twice continuously dif- 
ferentiate vector function with real- valued entries. The Pauli operator is 
given by the differential expression 



(3.12) Z = Q(a)®l + [. ^ "^.i + a^ 

+ 0-2,3 —0^3 



+ V® 1, 
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where 1 is the identity on C 2 , V = V(x) is the multiplication by a real- 
valued scalar potential and 

a i,k = ^ -x— , kj= 1,2,3. 

OXfc OXj 

Let B(x) be the length of the vector S(x) = curl a(x). Then the two 



eigenvalues of the perturbation of the term Q(a) <E> 1 in ( |3.12| ) at some 
point x G IR 3 are given by 

V(x) ±B(x). 

If V, B G U+! (IR 3 ) for some y > 3/2, then Theorem implies 



(3.13) trZl<L 



Y ^ T cl 

Y,3 



Q{(V + B)^ + (V-B)^}dx) 
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